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ABSTRACT 

We have discovered an optically rich galaxy cluster at z = 1.7089 with star formation occurring in 
close proximity to the central galaxy. The system, SpARCS104922.6+564032.5, was detected within 
the Spitzer Adaptation of the red-sequence Cluster Survey, (SpARCS), and confirmed through Keck- 
MOSFIRE spectroscopy. The rest-frame optical richness of A®gai(500kpc) = 30±8 implies a total 
halo mass, within 500kpc, of ^ 3.8+1.2x10®® Mq, comparable to other clusters at or above this 
redshift. There is a wealth of ancillary data available, including Canada-France-Hawaii Telescope 
optical, UKIRT-K, Spitzer-IRAC/MIPS, and Herschel-SPIRE. This work adds submillimeter imaging 
with the SCUBA2 camera on the James Clerk Maxwell Telescope and near-infrared imaging with the 
Hubble Space Telescope (HST). The mid/far-infrared (M/FIR) data detect an Ultra-luminous Infrared 
Galaxy spatially coincident with the central galaxy, with Lir = 6.2+0.9x10®^ Lq. The detection of 
polycyclic aromatic hydrocarbons (PAHs) at z = 1.7 in a Spitzer-IRS spectrum of the source implies 
the FIR luminosity is dominated by star formation (an Active Galactic Nucleus contribution of 20%) 
with a rate of ~86O+13OM0 yr“®. The optical source corresponding to the IR emission is likely a 
chain of of > 10 individual clumps arranged as “beads on a string” over a linear scale of 66 kpc. 

Its morphology and proximity to the Brightest Gluster Galaxy imply a gas-rich interaction at the 
center of the cluster triggered the star formation. This system indicates that wet mergers may be an 
important process in forming the stellar mass of BCGs at early times. 

Subject headings: galaxies: interactions, galaxies, starburst, galaxies: evolution, galaxies: clusters: 
general, infrared: galaxies 
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1. INTRODUCTION 

The most massive galaxies in the local universe reside 
at the centers of galaxy clusters. Often called Bright¬ 
est Gluster Galaxies (BGGs), they exhibit, as a class, 
highly uniform properties such as luminosity, metallicity 
and surface brightness profiles, and are distinct from the 
general galaxy population (Oemler 1976; Dressier 1978; 
Schombert 1986). These facts suggest formation pro¬ 
cesses which are themselves distinct from those which 
dominate in massive galaxies elsewhere in the cluster or 
in the field. The mass growth of BGGs is likely linked 
to the overall physics of hierarchical structure formation 
on galaxy cluster scales, including the fundamental pro¬ 
cesses of gas cooling, star formation, energy feedback and 
galaxy mergers at the centers of giant dark matter halos. 

Within the current suite of Semi-Analytic Models 
(SAMs) BGG evolution follows a relatively simple pre¬ 
scription (Kauffmann et al. 2000; Springel et al. 2001, 
2005; De Lucia et al. 2004, 2006; Croton et al. 2006). 
At very early times (z > 5) cooling flows dominate 
the stellar mass assembly at the centers of clusters, but 
these are rapidly suppressed through active galactic nu¬ 
clei (AGN) feedback (Bower et al. 2006; De Lucia & 
Blaizot 2007; Guo 2010). BCG growth then continues 
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over time through the accretion of smaller systems, with 
little additional star formation. In this picture most of 
the stellar mass is formed at very high redshifts {z ^ 5) 
within multiple galaxies through low levels of star forma¬ 
tion, and the BCG slowly acquires its identity through 
the conglomeration of previously assembled pieces. This 
is supported by direct measurements of the mass growth 
of BCGs below z ~1 (Lidman et al. 2012, 2013). 

Still, there are observational examples of moderate red- 
shift BCGs {z < 0.6) that present signs of high star for¬ 
mation rates (SFRs) (lO-lOOOM0yr“^) and/or contain 
large amounts of molecular gas (Allen et al. 1992; Mc¬ 
Namara et al. 1992; Jonhstone et al. 1998; Edge 2001; 
Wihnan et al. 2006; Edwards et al. 2007; O’Dea et al. 
2010; Hicks et al. 2010; Donahue et al. 2011; Rawle et al. 
2012; McDonald et al. 2012). Such active BCGs are rare 
and have primarily been seen in X-ray-selected cool-core 
clusters. In many cases the SERs are correlated with 
the gas cooling time. The paucity of these systems at 
low/moderate redshifts indicates that explosive star for¬ 
mation, fed by the cooling of halo gas, is not an impor¬ 
tant contributor to BGG mass growth, relative to the 
process of dry merging. 

While it is established that a large fraction (50%; Lid- 
man et al. 2013) of the stellar mass of BCGs is already 
in place at 2; ~ 1, it is not yet determined if this mass 
was assembled through processes similar to those at lower 
redshift (i.e., dry merging) or if a different and perhaps 
more extreme phase of growth, such as that seen in some 
of the cool core clusters, becomes more common at ear¬ 
lier epochs. There are now several detections of clusters 
and massive groups at redshifts beyond z 1 (Andreon 
et al. 2009; Muzzin et al. 2009; Wilson et al. 2009; Stan¬ 
ford et al. 2012; Zeimann et al. 2012; Gobat et al. 2013; 
Muzzin et al. 2013; Yuan et al. 2014; Bayliss et al. 2014; 
Tozzi et al. 2015) and these systems offer the opportunity 
to study the build-up of the central stellar mass in de¬ 
tail. Indeed many of these systems exhibit an increase in 
the star formation activity within several hundred kpc of 
their centers, compared to lower redshifts and, in some 
cases, the field (Tran et al. 2010; Brodwin et al. 2013; 
Alberts et al. 2014; Santos et al. 2015). While studies 
of their BCGs, in particular, have not been undertaken, 
these results indicate that the cores of galaxy clusters 
may be active places beyond z ~ 1.5. 

Here we report the discovery of an optically rich cluster 
at z =1.7089, that we label SpARCS104922.6-k564032.5 
(hererafter SpARGS1049-|-56). This system is distin¬ 
guished by its extremely active BGG which is observed 
in an important phase of evolution. We appear to 
be witnessing the stellar mass build-up of the central 
BGG through a wet galaxy merger, a process which has, 
heretofore, not been thought to play an important role in 
BGG growth. In this paper we outline the discovery and 
characterization of this interesting system. In Section 2 
we describe the data and observations; in Section 3 we 
present the analysis of the cluster and BCG properties; 
and in Section 4 we discuss the implication of these first 
results. We use standard cosmology throughout (Do = 1, 
Da = 0.7, Ho = 75 km s“^Mpc“^), and AB magnitudes. 

2. DATA AND OBSERVATIONS 

The data used in this work were obtained with a broad 
range of facilities. Some were retrieved from public 


archives while others were newly acquired by our team. 
Here we outline three categories of observations. In Sec¬ 
tion 2.1 we list the wide-held imaging associated with 
large surveys of the Spitzer- Wide-area InfraRed Extra- 
galactic (SWIRE) Lockman held, where the cluster is lo¬ 
cated. In section 2.2 we describe follow-up imaging with 
SCUBA2 on the James Clerk Maxwell Telescope (JCMT) 
and WEC3 on the Hubble Space Telescope (HST) taken 
by our team subsequent to the cluster discovery (de¬ 
scribed in Section 3.1). These data are focused on the 
location of the cluster and do not cover the entire Lock- 
man held. In Section 2.3 we present the spectroscopic 
observations obtained by others before the discovery of 
the cluster and by ourselves following the discovery with 
Keck-MOSFIRE. 

2.1. Existing Wide-field Imaging of the 
Lockman-SWIRE Field 

2.1.1. The Spitzer Adaptation of the Red-sequence 
Cluster Survey: 

The cluster was detected within the Spitzer Adaptation 
of the red-sequence Cluster Survey (SpARCS) coverage 
of the Lockman Hole. SpARCS is a near-infrared survey 
for galaxy clusters beyond z ~ 1 (Muzzin et al. 2009; Wil¬ 
son et al. 2009; Demarco et al. 2010). It combines deep 
2'—band imaging (5cr z' = 24.2 AB with the Canada- 
France-Hawaii Telescope; CFHT) with the SWIRE Sur¬ 
vey (Lonsdale et al. 2003) at 3.6^m, to identify over¬ 
densities of red-sequence galaxies at high redshift. A 
thorough discussion of the z' data analysis, z' — 3.6/rm 
cluster hnding algorithm and the SpARCS cluster cata¬ 
log will be presented in a future paper by A. Muzzin et 
al. (2015, in preparation), though we outline the method 
in Section 3.1. 

2.1.2. The Spitzer Space Telescope: 

The public SWIRE Legacy Survey consists of imaging 
at 3.6, 4.5, 5.8 and 8.0/im with IRAC and 24, 70, 160/xm 
with MIPS (Lonsdale et al. 2003) SpARCS1049-l-56 

lies near the edge of the survey and is not covered by the 
SWIRE ancillary optical data for this field. Deeper (~ 
2/rJy) 3.6/4.5/rm imaging of this area has been obtained 
during the warm-Spitzer mission by the Spitzer Extra- 
galactic Representative Volume Survey (Mauduit et al. 
2012), but it is not used here. 

2.1.3. The Herschel Multi-tier Extragalactic Survey: 

The Lockman-SWIRE field is also part of the Herschel 
Mid-Infrared Extragalactic Survey (HerMES; Oliver et 
al. 2012). The primary data for this survey are SPIRE 
250/350/500/rm imaging (Griffin et al. 20120), with ad¬ 
ditional PACS data taken in parallel mode. PAGS data 
products have not yet been released by the HerMES team 
and therefore we downloaded the original data from the 
Herschel Science Archive. We produced maps at 100 and 
160/xm following the technique developed in Ibar et al. 
( 2010 ). 

2.1.4. The UK Infrared Deep Sky Survey: 

http : //swire.ipac.caltech.edu/swire/astronomers/publications/ 
SWIRE2_doc_083105.pdt 
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This field was imaged in the near-infrared as part of the 
UK Infrared Deep Sky Survey (UKIDSS). The UKIDSS 
project is defined in Lawrence et al (2007). UKIDSS uses 
the UKIRT Wide Field Camera (Casali et al, 2007). The 
photometric system is described in Hewett et al. (2006), 
and the calibration is described in Hodgkin et al. (2009). 
The pipeline processing and science archive are described 
in M. Irwin et al (in preparation) and Hambly et al 
(2008). The Lockman Hole is contained within the Deep 
Extragalactic Survey and reaches a depth of K=21. As it 
is located near the edge of the survey, SpARCSI049-|-56 
is only covered by the K-band. 

2.1.5. New CFHT Imaging of the Lockman Field 

As part of SpARCS follow-up imaging, 35 square de¬ 
grees overlapping the Spitzer Adaptation of the Red- 
sequence Cluster Survey Northern Fields were observed 
in the ugrz' bands using MegaCam on the 3.6 m CFHT, 
including the Lockman field (A. Tudorica et al. in prepa¬ 
ration). Exposure times were 4201, 2291, 3007, and 6001 
s for the urgz' filters respectively, and seeing averaged 
0.8 arcsec. The data were reduced using the ELIXIR 
and THELI pipelines. The absolute photometric cali¬ 
bration is based on the SDSS DRIO, which was used 
to determine the field-by-field zeropoint values. The 5a 
limiting magnitudes are u = 25.Qmag,g = 25.3mag,r = 
25.0mag, z' = 2AAmag. Eor this work these data are 
used for display purposes only, in Figure 2. 

2.2. New Imaging of SpARCS1049+56 
2.2.1. JCMT SCUBA2 450/850pLm 

We obtained 3 hours of integration in Band-1/2 
weather with the new SCUBA2 bolometer array on the 
JCMT. SCUBA2 simultaneously observes at 450/im and 
850/rm with 8" and 15" FWHM beams respectively. We 
employed the simple Daisy scan strategy which provides 
uniform sensitivity over the center ^3' area. The maps 
were produced using the standard SCUBA2 data reduc¬ 
tion pipeline (DRP) (Jenness et al. 2011; Geach et al. 
2013; Chapin et al. 2013) and beam-filtered to maximize 
sensitivity to point-sources. In the central region of each 
map (located at the cluster core) we reach an rms of ^ 8 
mjy beam“^ at 450/im and ^ ImJy beam“^ at 850/im. 

2.2.2. WFC8 Imaging 

We obtained infrared imaging of the cluster center with 
the F160W and F105W filters in Cycle 22 These 
observations were obtained over nine separate visits, with 
total exposure times of 9237s and 8543s in F160W and 
F105W respectively. 

The data were processed using standard HST data re¬ 
duction techniques, with some minor modifications. To 
combine and re-sample all images, we use the standard 
Space Telescope Science Institute (STSci) AstroDrizzle 
software from the DrizzlePac package , using the “ime- 
dian” combine type with 0.09 pixels. The DrizzlePac 

Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained [from the Data Archive] at the Space 
Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract 
NAS 5-26555. These observations are associated with program # 
GO-13677 and # GO-13747 

http://drizzlepac.stsci.edu/ 


software accounts for all distortion corrections, and also 
accounts for other instrument/detector specific calibra¬ 
tions, much of which is contained in the image data qual¬ 
ity array. Cosmic rays are flagged and corrected in the 
up-the-ramp sampling, and our 3-dither pattern allows 
unfiagged hot pixels or uncorrected cosmic rays to be 
found and masked by AstroDrizzle. 

To align the images, we detect objects in each image 
with SEP and then apply standard point-set regis¬ 
tration techniques to shift and rotate the images. We 
also use SEP to perform spatially variable background 
subtraction, masking the bright objects in the image to 
avoid bias in this background estimate. 

2.3. Spectroscopy 

2.3.1. Keck-MOSFIRE Spectroscopy 

H-band spectroscopy was obtained in 2013/14 using 
the near-infrared multi-object spectrograph MOSFIRE 
on the W.M. Keck I Telescope (McLean et al. 2012). The 
MOSEIRE H-band filter is centered at 1.632//m, and has 
a spectral resolution R = 3660. This allows the resolved 
observation of Ha and the [Nn] doublet, AA6583, 6548 
at the redshift of SpARCS1049-|-56. Pour separate H- 
band masks were obtained with ^ 40 slits per mask. 
During observations, each individual exposure was 120s 
in length, with total integration times 45-60 minutes per 
mask. All observations were taken utilizing an ABBA 
dither pattern with a nod amplitude of 1.5 and standard 
0/7 slit width. 

We utilized the MOSFIRE data reduction pipeline 
(DRP v20140610) for the first step in data reduction. 
The DRP produces a two-dimensional spectrum and a 
corresponding sigma image from the raw flat and data 
frames for each object on the slit mask. All spectrum 
and sigma images are a sky-subtracted, wavelength cali¬ 
brated, then registered to produce a stack of the A and B 
dither positions. All raw images were visually inspected 
for integrity before ingested into the DRP; any anoma¬ 
lous frames were thrown out. One step in the DRP in¬ 
volves interactively fitting the wavelength solution based 
on the sky lines during which we attempt to fit the wave¬ 
length solution to have a maximum residual difference of 
0.2 pixels for any one skyline. 

We then extracted one-dimensional spectra from the 
two-dimensional spectra and sigma images following the 
extraction method of Horne (1986). Eull details will be 
present in a future paper by DeGroot et al. (2015, in 
preparation) but we give a brief outline here. Typical 
extraction routines simply sum the the sky-subtracted 
image data from the two-dimensional spectrum over a 
range that encloses the object’s spectrum in the spatial 
direction. Our extraction applies a non-uniform pixel 
weighting, nominally gaussian in shape, to the extrac¬ 
tion sum which increases the signal-to-noise of the result¬ 
ing one-dimensional spectra. Redshifts are determined 
from developed line fitting routines that simultaneously 
fit Ha and [Nn] doublet emission lines as well as a linear 
continuum component via a Markov Chain Monte Carlo 
(MCMC) approach based on bayesian inference. Red- 
shift errors are determined from the addition in quadra¬ 
ture of the random error from the MCMC fit and the 

https://github.com/kbarbary/sep 
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systematic error from the interactive wavelength fitting 
step of the DRP. 

2.3.2. Spitzer IRS Spectroscopy 

The position of the galaxy that we identify as the cen¬ 
tral BCG (Section 4.2) was serendipitously observed as 
part of an earlier Spitzer-IRS survey of Farrah et al. 
(2008) (see Section 3.3 for further discussion of the identi¬ 
fication of the BCG with this observation). This program 
targeted infrared-bright systems (S 24 ;im > SOOmJy) with 
infrared colors indicative of high-redshift galaxies. The 
data acquisition, reduction and analysis are described in 
Farrah et ah, and we use these measurements in Section 
3.3. We further downloaded the spectrum from the Cor¬ 
nell Atlas of Spitzer IRS Spectra (CASSIS) which has 
developed its own analysis pipeline (Lebouteiller et al 
2011 ). 

3. ANALYSIS 

3.1. Diseovery and Confirmation of the z = 1.7 
SpARCS1049+ 56 cluster 

The northern/equatorial SpARCS SWIRE images 
cover 9.8, 4.5, 11.6 and 9.4 deg^ (for a total of ^ 35 deg^) 
over the ELAIS-Nl, ELA1S-N2, Lockman and XMM-LSS 
multi-wavelength survey fields respectively (Muzzin et al. 
2009). SpARCS1049+56 is located within the Lockman 
field (the largest single field) and is the second richest 
cluster detection above z = 1.7 within the 35 square de¬ 
grees. Unfortunately, it lies close to a chip gap on the 
z' image and the outskirts of the northern half of the 
cluster lack coverage in this band. 

The Keck-MOSEIRE spectroscopy confirms the pres¬ 
ence of a galaxy over-density at z = 1.7. An average 
around the z ^ 1.7 peak produces a cluster redshift of 
z = 1.7089. In Figure 1 we show the distribution of 70 
redshifts measured from four MOSFIRE masks of 146 in¬ 
dependent slits, over a field-of-view of ^ 6x3 arcmin, and 
in Figure 2 we show the physical locations of the galax¬ 
ies with spectra. The distribution of confirmed cluster 
members (see discussion below) is spatially extended over 
several Mpcs, but this is, in part, a consequence of our 
observing constraints rather than a reflection of the true 
spatial distribution. The MOSFIRE design prohibits 
complete sampling of the core, due to a required min¬ 
imum separation between slits, and thus the confirmed 
members are spatially spread out (Figure 2). Moreover, 
redshifts were obtained through the detection of Ha and 
this could impose a radial bias in the redshift success, 
since completely passive galaxies (which may dominate 
within the core) will not show emission lines. 

We estimate the cluster richness from the 3.6/im (rest- 
frame ^J) data to be iVgai =30 ± 8 . We use the fol¬ 
lowing definition of AgaP the number of background- 
subtracted galaxies within a 500kpc radius of the cluster 
center brighter than M*+\ based on a passive evolving 
luminosity function with zformation = 10 (Muzzin et al. 
2007). Applying the conversion of Andreon & Congdon 
(2014b) this richness yields a mass within 500 kpc of 
(3.8±1.2)Xwhere the uncertainty does not in¬ 
clude redshift evolution in this relation. We note that 
in a richness-limited sample, such as this, there may be 
a bias toward systems with a high richness-to-mass re¬ 
lation. Confirmation of this mass will require a gravita- 



Fig. 1.— The distribution of redshifts from the Keck-MOSFIRE 
spectroscopy. The inset shows a zoom in on the redshift peak of 
z = 1.7089. A velocity cut of ±1500 km/s of the mean cluster 
redshift results in 27 members within the 6x3 arcmin field of view 
of MOSFIRE. 


tional leasing measurement or detection of the intraclus¬ 
ter medium (ICM). Currently there is no X-ray coverage 
of this object, and it is not expected (given its estimated 
mass) to appear as an SZ decrement in the Planck all-sky 
map (Planck Collaboration et al. 2014). 

We classify as cluster members galaxies which lie 
within ±1500 km/s of the mean cluster redshift. This 
provides a sample of 27 cluster members within a cluster¬ 
centric radius of 1.8 Mpc. A standard biweight-mean 
yields a velocity dispersion of 430tioo km/s. We use 
the model calibrated relation between cluster viral mass 
and velocity dispersion to convert this to a virial mass of 
8±2xl0^^ Mq. However, Saro et a. (2013) have found 
that calculating the velocity dispersion using at least the 
brightest 30 red sequence galaxies within the virial radius 
is required to obtain the virial mass of a cluster to within 
40%, whereas we have used 27 Ha emitters over a wider 
area. This region of the spectrum is also plagued by sky 
lines and does not allow us to fully, or symmetrically 
sample the redshift distribution. We therefore attach an 
uncertainty of at least 40% to this mass, or 8 ± 3x10^^ 
Mq. 

Although the mass is uncertain, we can add this sys¬ 
tem to the growing list of clusters and massive groups 
known around or above z = 1.7. The most massive 
published clusters are the Stanford et al. (2012) sys¬ 
tem at z = 1.75 (IDCSJ1426±3508) with seven spec¬ 
troscopically confirmed members and X-ray and SZ de¬ 
tections that indicate a mass estimate of M 200 '^5x10^'’’ 
Mq (Brodwin et al. 2012), and the Tozzi et al. (2015) 
cluster with M 500 =3.2f_Q'^xlO^'^MQ. JKCS041 (An¬ 
dreon et al. 2009, 2014a; Newman et al. 2014) lies at 
z =1.803, with a mass constrained by optical richness 
and X-ray properties to be 1 . 6 xlO^^M 0 . The Zeimann 
et al. ( 2012 ) system is at a slightly higher redshift of 
z = 1.89, but with no mass estimate reported. Finally, 
CL1449±0586 (Gobat et al. 2013) is confirmed at z = 2.0 
with 26 spectroscopic members. Its X-ray flux indicates 
a mass of 5 x10^^Mq. With a mass of 0.8-3 .8 xlO^^ Mq 
SpARCS1049 is similar to these objects. 
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Fig. 2. — Left: The X— band image (greyscale) of the ^^6x6 arcmin region around SpARCS1049H-56. Circles denote galaxies with spec¬ 
troscopy from Keck-MOSFIRE. Large red circles indicate confirmed cluster members; smaller dark blue circles show foreground/background 
objects; and smaller cyan circles correspond to galaxies for which no redshift could be obtained. This latter group contains fore¬ 
ground/background systems as well as cluster members with no strong emission lines. The central green circle shows the location of 
the BCG, although the BCG itself is not easily visible in the this image (as it is very faint and red). Right: The ur3.Q{im image of the 
cluster region with the location of the BCG again indicated. 
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3.2. The Galaxy Population and the Brightest Cluster 

Galaxy 

Figure 3 presents the color-magnitude diagram for 
SpARCS1049-f56 using HST-WFC3 NIR filters which 
bracket the 4000A-break in the rest-frame. We show all 
objects detected in both hlters within the central 40" (~ 
320 kpc) radius of the HST image. Orange stars show 
galaxies confirmed as foreground or background systems, 
with the largest star corresponding to a low-redshift 
line-of-sight companion to the BCG (J.S. Fluang, pri¬ 
vate communication). Solid blue circles indicate cluster 
members confirmed through emission-line spectroscopy. 
Large open red circles correspond to galaxies for which a 
spectrum was obtained, but no redshift was determined, 
due to a lack of features. These include fore/background 
systems and cluster members with no strong emission 
lines in the iJ-band, and emission-line cluster galaxies 
whose Ha line is masked by a sky-line. 

An over-density of red-sequence galaxies is observed; 
this is expected since that is required for the initial 
cluster detection. We show the location of the red- 
sequence (black line) predicted for a single stellar pop¬ 
ulation (SSP) formed at Zmathrmform = 5 (Bruzual 
& Chariot 2003). It is interesting to note that all of 
the galaxies with colors consistent with the red sequence 
do not show evidence of strong line emission, implying 
the Red Sequence is already old and passive in this sys¬ 
tem. Nor are the identified cluster members very blue, 
rather they clump just below the red-sequence in the ap¬ 
proximate region of the “green-valley”. This might be 
explained by the two primary biases of the MOSFIRE 
spectroscopic observations. Targets were prioritized by 
z' - 3.6^m color, to best select cluster members. This 
biases the spectroscopy toward red galaxies. However, 
redshifts required the detection of the Ha emission line, 
thereby requiring star formation. This combined selec¬ 
tion could tilt the confirmed members toward the reddest 
emission-line galaxies, or those in the green valley. Still, 
the presence of these systems may indicate that a sig¬ 
nificant fraction of the galaxies within SpARCS1049 are 
transitioning from a star-forming episode to passive red- 
sequence galaxies. 

The BCG stands out clearly on the Red Sequence and 
is labeled in Figure 3 with the next brightest system al¬ 
most a full magnitude fainter. Ha was not observed for 
this galaxy, indicating it has a passively evolved stel¬ 
lar population (and see Section 3.3) with no measurable 
trace of star formation. However, the IRAC colors in¬ 
dicate, through the presence of the stellar photospheric 
feature at rest-frame 1.6/xm (the stellar bump), a photo¬ 
metric redshift of z = 1.7. 

3.3. Ultraluminous Infrared Emission Coincident with 

the Central Galaxy 

The center of the cluster is coincident with a previ¬ 
ously known Ultraluminous Infrared Galaxy at z = 1.7 
in the Farrah et al.(2012) IRS study of extremely IR- 
bright high-redshift galaxies in SWIRE. The sample was 
defined by IRAC colors, to pre-select galaxies at z ~ 1.7, 
and bright MIPS-24/rm (> 500/xJy) emission. The target 
of the Spitzer spectroscopy was the BCG of the cluster 
and strong polycyclic aromatic hydrocarbon (PAH) fea¬ 
tures were detected, confirming a redshift of z = 1.7. 



Fig. 3. — The HST NIR color-magnitude diagram, for all ob¬ 
jects within 40” of the cluster center. For clarity we do not include 
photometric uncertainties, but show the ~5 (t color-limit set by the 
F105W limit (diagonal dotted line); the same limit for F160W is 
~ 25.5 mag. Small open circles correspond to all objects with 
photometric measurements in the inner 40” radius of the HST 
image. Solid orange stars indicate confirmed non-cluster mem¬ 
bers (generally foreground), with the largest star corresponding to 
the bright blue line-of-sight companion to the BCG (J.S. Huang, 
private communication). Larger red circles highlight galaxies for 
which a spectrum was obtained, but no redshift was determined 
because of a lack of spectral features. These are a combination of 
fore/background galaxies, and galaxies within the cluster with no 
strong emission lines (passive systems). Finally, the large, solid 
blue circles indicate cluster members confirmed through emission¬ 
line redshifts. We also include a model of the Red Sequence for 
a formation redshift of Zform = 5. The identified BCG (labeled) 
stands out as the brightest galaxy on the red-sequence. 

Two other bright galaxies can be seen within the MIPS 
24/rm beam, to the SE. Both of these galaxies have been 
spectroscopically confirmed to be foreground objects at 
z = 0.64 (J.S. Huang, private communication, the Multi- 
Mirror Telescope ) and z = 1.3 (Keck-MOSFIRE), for 
the blue and red galaxy respectively. The detection of 
PAH lines at z = 1.7 therefore removes these two galaxies 
as candidates for the MIPS emission. Thus, the ground 
based-imaging leaves only the BGG as the candidate ID; 
although some diffuse emission is visible, no other source 
is seen within the 95% confidence region. 

The deep HST imaging reveals a more complicated 
situation, however. In Figure 4 we show the two-color 
composite image of the central region using the F105W 
and F160W filters (rest-frame ^ u' and r', respectively). 
Here we see a long chain of small objects embedded 
within a diffuse filament or tidal-tail. The MIPS centroid 
is located directly on top of one of these objects (in the 
approximate center of the chain) and thus it seems likely 
that the MIPS emission is coming not from the BGG it¬ 
self, but from an object, or set of objects, within ~16kpc 
of the BCG. The details of this structure are discussed in 
the following section, but we note its morphology and the 
proximity to the BCG (possibly within its stellar enve¬ 
lope) suggest that the strong infrared emission is driven 
by an interaction with the central galaxy of the cluster. 

3.3.1. Constraining the Infrared Spectral Energy 
Distribution 

With such extensive multi-wavelength coverage of 
infrared-bright object, we are in a position to constrain 
its spectral energy distribution (SED) and derive sim- 



7 


1/3 


I 


C 

o 


10 5 0 -5 -10 

Arc Seconds 



Fig. 4.— Top: the two-color (F160W, F105W) HST image of 
the central region of SpARCS1049, showing the complex morphol¬ 
ogy of the core. Diffuse emission in a tidal-tail shaped structure is 
visible, and embedded along its length are multiple clumps. The 
cross marks the centroid of the 24/i.m emission. Bottom: The 
3.6/i.m image of the same region (greyscale) with MIPS flux con¬ 
tours overlaid (orange) to show the location of the 24/im centroid 
and the mild extension of the source in the N/S direction. The two 
purple crosses correspond to the 850/i.m (large) and 450/im (small) 
positions, where the size of the cross denotes 1/2 the FWHM of 
the beams. 

pie parameters from it. Flux measurements in the opti- 
cal/NIR/MIR regimes are straightforward, however the 
large beam sizes of the Herschel instruments make confu¬ 
sion difficult to correct for in the crowded central regions 
of the cluster. We therefore explain our method in detail 
here. 

The SPIRE data are publicly available and while there 
is emission coincident with BCG in the HerMES 250^m 
catalog (Wang et al. 2013), the increasing beam-size at 
350 and 500/j,m results in a highly confused image of 
the cluster core. However, utilizing positional priors 
from less-confused wavelengths can reduce the amount 
of source blending in far-infrared photometry. In par¬ 
ticular, Viero et al. (2013) developed a novel algorithm 
called SIMSTACK to stack SPIRE maps on prior posi¬ 
tions from a user-defined catalog. SIMSTACK simulta¬ 
neously stacks lists of potentially clustered populations. 


which effectively allows for a deconvolution of fluxes from 
multiple sources within each SPIRE beam. We exploit 
the S'pjfzer-MIPS 24/rm catalog from SWIRE for our po¬ 
sitional priors. 

As we are only interested in the flux of the infrared 
source , we alter the code slightly to suit our needs. We 
collect all 24/rm sources within 90 arcsec of the BCG, 
which corresponds to ~50, 30, and 15 independent beams 
at 250, 350, and 500/im, respectively. We separate each 
source into individual lists, and simultaneously fit each 
list for the best value, yielding a flux of 46 ±6.2, 47±6.7, 
and 21 ±9.4 for ascending SPIRE wavelengths, including 
the nominal confusion noise (Nguyen et al. 2010). We 
compare these fluxes to those derived from the HerMES 
XID catalog that uses the 250/im source positions (which 
have a FWHM beam ^ 3 times larger than 24^m) as 
priors for 350 and 500/j,m. While the individual SPIRE 
fluxes can differ by 2 — 3 (t between the two methods, 
the resulting best fit SED (see Figure 5 and below), and 
subsequent infrared luminosity, are robust, differing only 
by 1.3%. 

At 850^m we detect a point source at an offset of 4" 
from the BCG location with a significance of A.2a. A 
450^m peak is detected at the same location at signifi¬ 
cance of 4.0(7 (Figure 5). Both of these offsets are within 
the expected positional error radius for SCUBA2, given 
the signal-to-noise of the detection, and these measure¬ 
ments are included in Figure 5. No useful constraints 
are obtained from the MIPS 70 and 170/im imaging of 
SWIRE. 

In Figure 5 we show the complete SED for the far- 
infrared source using all available data from 24 to 850^m 
(Table 2), including the IRS spectrum in cyan. The re¬ 
duced IRS spectrum was obtained through the CASSIS. 
We have ht two sets of template SEDs to the data: we 
use the the composite SED of Kirkpatrick et al. (2012) 
of star forming ULIRGs at ^ ~ 2, and the library of 
templates from Chary & Elbaz (2001). We fit only to 
the photometry, for consistency, but include a discussion 
of the information contained within the IRS spectrum 
below. 

The Kirkpatrick SED is derived from infrared- 
luminous objects of similar luminosities to this source, 
at a similar redshift of z ~ 2. As there is only one spec¬ 
trum, we fit for the amplitude and then integrate under 
the infrared portion (8-1000/im) of the SED to obtain 
the total infrared luminosity. Uncertainties on the lumi¬ 
nosity are produced by htting the extremes of the uncer¬ 
tainties/scatter on the composite template. 

The Chary & Elbaz templates, on the other hand, are 
derived from local galaxy templates and thus may not be 
appropriate for galaxies at z 2. Indeed, many studies 
have shown that using the CEOl SEDs to infer the Lir 
from a single 24^m measurement over estimates the lu¬ 
minosity beyond z = 1.5 (Murphy et al. 2009; Nordon et 
al. 2012; Rodrghiego et al. 2010). As explored by (Elbaz 
et al. 2011), however, this discrepancy is reduced when 
one considers the full range of CEOl templates, rather 
than pre-selecting a single template based on the ob¬ 
served 24^m luminosity, and when fitting broader photo¬ 
metric coverage. We therefore adopt the following fitting 
method for these templates. We fit each available tem¬ 
plate (105 total variations) to the data (Pope et a. 2008). 
We calculate uncertainties on the infrared luminosities 





using a Bayesian approach. We calculate the normalized 
chi-square probability distribution as a function of Chary 
& Elbaz template and amplitude, ranging from 0-50 in 
steps of 0.05. Each value in this grid corresponds to an 
infrared luminosity, which can then be weighted by the 
normalized chi-square probability. The final uncertainty 
is given by the standard deviation of the weighted mean 
of the infrared luminosity: ctlir = \J{LIR?) — (LIR)'^. 

In Table 2 we summarize the results of these fits. For¬ 
mally, the best representation of the data is the CEOl 
SED This fit produces Lir = 6.2±0.9xl0^^ Lq, with a 
xt = 0.4. We note the uncertainty on this value is small 
because only a small set of CEOl SEDs accurately repre¬ 
sent the data. The Kirkpatrick model yields Lir=4.1± 
1.5x10^^ Lq, with xi = 4.0. As can be seen in Fig¬ 
ure 5, the Kirkpatrick model over-predicts the amount 
of long-wavelength emission in the Rayleigh-Jeans tail. 
However, the integrated luminosities of the two fits are 
consistent with each other at a level of ^ Ict. 

The LIR emission can be converted to a SFR by ap¬ 
plying the standard Kennicutt (1998) law, under the as¬ 
sumption that the dust is heated entirely by the stellar 
population with a Salpeter IMF, and these values are 
also included in Table 2. The CEOl fit yields a SFR of 
1070±160 M 0 yr“^, while the Kirkpatrick fit produces 
690±260 M 0 yr“^. These estimates are in agreement 
with the upper-limit on the SFR placed by Farrah et 
al. (2008) of < 1200 M 0 yr“^ from the non-detection of 
the 6.2 fim line. 

If an AGN is contributing to the IR luminosity this 
estimate must be reduced, however the presence of the 
PAH features in the spectrum, and the turn-over of the 
IRAC fluxes (implying a lack of very hot dust), suggest 
that an AGN does not dominate the luminosity. To in¬ 
vestigate this we decompose the IRS spectrum into a star 
forming and AGN components following the method of 
Pope et a. (2008). We find an AGN contribution to the 
MIR of 50%, corresponding to an AGN contribution to 
the FIR luminosity of ^^20%. We therefore apply this 
downward correction to the SFR. This yields, for the 
best-fit SED, a SFR of 860±130 Mgyr-b 


3.4. HST-revealed Beads on a String 

Prior to the HST NIR data, the BGG was the best 
candidate for the 24^m emission as it appeared to be the 
only viable option within the search radius of the source; 
however, the high resolution of the HST image has re¬ 
vealed a very complex structure within the core of the 
cluster, and has detected numerous additional objects 
within the search area. A color zoom-in of the central 
region is shown in Figure 4. Here we see a chain of small 
clumps arranged as “beads on a string” and covering a 
linear area of roughly 66 kpc. The clumps are embed¬ 
ded within diffuse emission which has the morphology of 
a single large tidal tail. This tail appears to begin 1" 
east (8 kpc) of the BCG, within its stellar envelope, and 
forms a distinct backwards J-like structure, ending to the 
south-west of the BCG. 

The morphology of this structure is reminiscent of in¬ 
teracting galaxies. In particular, it is similar to the long 
single tidal arms of the “shrimp” galaxies at z ^ 1.4 
classified by Elmegreen et al. (2007) in the GEMS and 


GOODS fields. These systems appear to be higher red- 
shift analogs of merging systems such as M51 or the An¬ 
tennae. They consist of a single large tail, defined by 
diffuse emission and composed of several distinct and 
regularly spaced clumps. The scale of the structure in 
SpARCS1049-|-56 is similar to these, with a linear extent 
of 66 kpc, a width of ^^5 kpc and clump spacing ranging 
from 1 to 2 kpc. The mips 24/rm emission peaks on one 
of the beads in the center of the tail and appears to be 
slightly extended along the direction of the tail (Figure 

4) 

4. DISCUSSION 

4.1. Dissipative Merger Induced Formation of the BCG 

The unique feature of this cluster is its extremely ac¬ 
tive and infrared-luminous central object. We conclude, 
based on the presence of PAH features and its NIR col¬ 
ors, that the IR emission is dominated by star formation 
rather than an AGN. The best SED fit yields Lir = 
6.2±0.9x 10^^ Lq and a corresponding AGN-corrected 
SFR of 860±130 MQ-yi~^. In this respect the system 
is similar to the Phoenix BCG (McDonald et al. 2012) 
within a z = 0.6 SPT cluster, but SpARCSI049-|-56 is 
seen at a significantly earlier cosmological epoch. In both 
of these systems a substantial amount of gas is being 
delivered to the cluster center to fuel this intense star- 
burst. In the case of the Phoenix the gas is deposited by 
a rapidly cooling intracluster medium (ICM) (McDonald 
et al. 2013). X-ray imaging of the Phoenix cluster indi¬ 
cates it has a cool core and deep optical and UV imaging 
and spectroscopy has revealed a network of large-scale 
gas filaments apparently streaming into the BCG (Mc¬ 
Donald et al. 2013, 2014). 

Although no X-ray information exists for 
SpARCS1049, the NIR (rest-frame UV and opti¬ 
cal) data indicate a very different scenario. Rather than 
a spider-web arrangement of filaments, the HST imaging 
shows a single tidal tail, with clumps along its length. 
While this does not rule out the deposition of gas by 
a cooling flow, the morphology is more consistent with 
starbursts induced by gas-rich mergers. 

Though the phenomenon of mergers producing “beads 
on a string” tidal tails is well-established to our knowl¬ 
edge (Elmegreen & Efremov 1996) SpARCS1049 is only 
the second example of such an event in the center of a 
galaxy cluster. Recently, a similar object was reported 
by Tremblay et al. (2014) in a lower-redshift {z = 0.3) 
strong-lensing cluster, SDSSJ1531-I-3414. This system 
is clearly associated with a major-merger consisting of 
the BCG and a massive companion. There are sim¬ 
ilarities between SpARCS1049 and SDSSJ1531-I-3414, 
namely the spacing of the beads ranges between 1 and 2.5 
kpc and the linear extent is several tens of kpc long (66 
kpc and 27 kpc respectively). However, there are some 
key differences as well. Firstly, a second massive galaxy, 
merging with the BCG is not obvious in SpARCS1049. 
The tidal tail appears to begin within the stellar envelope 
of the BCG with a smaller merging system that is spec¬ 
troscopically confirmed to be within the cluster sitting 
at its tip. The tail may be associated with this system 
which itself is interacting with the BCG, but does not 
have a comparable mass. Secondly, although the rest- 
frame filters are necessarily different the beads/clumps 
of SpARCS1049 do not appear uniform in color, unlike 



TABLE 1 

Luminosity and Star Formation Rate Estimates for the 
Infrared-Bright Central Galaxy System 
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Method/SED 

xi 

Air 

(1012 Lq) 

SFR 

{Mq yr-i) 

AGN removed SER^ 
(Mq yr-l) 

Chary Elbaz 

0.4 

6.2 ± 0.9 

1070 ± 160 

856 ± 128 

Kirkpatrick 

4.0 

4.1 ± 1.5 

690 ± 290 

552 ± 232 

6.2/im PAH^ 



< 1200 



^Assuming 20% of the FIR luminosity is produced by the AGN 
^Farrah et al. (2008) 
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Fig. 5.— The infrared spectral energy distribution from the combined MIPS/IRS/SPIRE/SCUBA2 flux measurements. The photometric 
points are shown in dark red and the IRS spectrum in cyan.The red arrows show the 3a limits of the PACS measurements. Overlaid (solid 
black line) is the infrared luminous template SED for z 2 ULIRGs developed by Kirkpatrick et al. (2012), with the scatter in this SED 
denoted by the grey shading. The blue curve corresponds to the best-fit SED from the (Chary & Elbaz 2001) template library, with the 
blue shading showing the family of templates that produce a <3a fit. Both SEDs are scaled to best-fit the photometric points only; the fits 
do not include the IRS spectrum which is simply over-plotted as observed. The Chary Sz Elbaz SED (without the optical points) yields the 
lowest provides an infrared luminosity of Ljr = 6.6x10^^ Lq. The green dotted line shows a modified black body (/3= 1.7) fit to 

the SPIRE and SCUBA2 photometry which provides yet another estimate of the infrared luminosity, that is independent of the assumed 
SED. The dust temperature is found to be T = 44 di 5 K and Ljr 4x10^^ Lq, in in good agreement with that of the full SED. The 
parameters of the template fits are listed in Table 1. 


TABLE 2 

Infrared Photometry Measurements 


24/im (mjy) 

250/im (mJy) 

350/im (mjy) 

450/im (mJy) 

500/im (mjy) 

850/im (mjy) 

0.645 ± 0.065 

46.2 ± 6.2 

46.8 ± 6.7 

33.2 ± 8.3 

20.7 ± 9.4 

4.5 ± 1.1 
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those of SDSSJ1531+3414. The clumps near the BCG- 
end of the tail are blue and similar in color to the merging 
pair at the tip, while the extended end of the tail has a 
color more consistent with the BCG itself. Given that 
the MIPS emission is coincident with this second half of 
the tail, this may be due to extreme dust extinction in 
this region. Finally, the scale of the star formation in 
SpARCS1049+56 is extremely different. The total SFR 
of SDSSJ1531+3414 (determined from Ha) is 5 MQyr—1, 
approximately two orders of magnitude smaller than we 
measure for SpARCS1049. 

Assuming that the stellar mass produced in the burst 
will eventually be incorporated into the BCG we can as¬ 
sess its importance through a rough estimate of the spe¬ 
cific SFR. For this we simply use the total SFR as seen in 
the IR and the existing approximate mass of the BCG. 
To estimate the stellar mass we employ the color depen¬ 
dent rest-frame K mass-to-light ratio (M/Lif) calculated 
by Bell et al. (2003), using the F160W-F105W color as 
an approximate rest-frame u-g. This yields stellar mass 
of 3±0.4xl0^^ Mq, where the error reflects a 10% uncer¬ 
tainty in the photometry (corrected to a Salpeter IMF 
(Raue & Meyer 2012)). Note that in the rest-frame K 
galaxies have a weak dependence of M/L on color: an 
uncertainty in color of ±0.5 would lead to an additional 
uncertainty of ±0.5x10^^ Mq. Thus, although it is un¬ 
dergoing a very rapid starburst of 860Mq yr“^ , the spe¬ 
cific SFR of the SpARCS1049±56 BCG is a moderate 
3-5 Gyr“^, in line with mass-extrapolated main sequence 
field galaxies at 2 ~ 2 (Daddi et al. 2007). This could 
imply that the mechanism of gas deposition is similar to 
that of the field at this redshift (major gas-rich mergers), 
or that galaxy self-regulation operates in the same way 
to limit the star formation efficiency, regardless of the fu¬ 
elling process. Nevertheless, the inferred SFR of 860Mq 
yr“i would double the stellar mass of the galaxy within 
roughly 400 million years therefore necessitating a rapid 
quenching mechanism. This is a similar timescale to that 
expected for a cooling flow of 100-500 Myr (Sarazin et al. 
1986), with subsequent shut-down by an emerging AGN, 
but is also consistent with the gas depletion timescales 
for ULIRG mergers (Kennicutt 1998). 

Regardless of the central fueling mechanism, this is 
an optically rich cluster caught in the process of form¬ 
ing a large amount of stellar mass within the halo of its 
BCG. SpARGS1049±56 therefore offers a new opportu¬ 
nity to study the interplay between the global cluster 
physics, be it the cooling of the gas halo or the inter¬ 
action between infalling galaxies and the ICM, and the 
evolution of the massive central galaxy within the system 
at a time when such processes may be more important to 
building clusters than today. The rapid growth of BGGs 
through large-scale cooling of the IGM is a key ingredient 
at early times in SAMs (e.g. De Lucia & Blaizot 2007), 
but the supporting observational evidence is slim, largely 
because of the dearth of high-redshift {z > 1.5) over¬ 
densities available for study. Evidence is mounting that 
growth of BGGS through dry mergers is an important 
process (e.g. Lidman et al. 2012), but these studies have 
been limited to z ^ 1. It may be that SpARGS1049±56 
is the first example of a new process, wet mergers of gas- 
rich galaxies, occurring in the core of a z ~1.7 cluster. 

5. CONCLUSIONS 


This paper presents the discovery and first analysis 
of a rich galaxy cluster at z = 1.708 with a remarkable 
star-forming central galaxy. While much additional work 
must be done to characterize the cluster itself we can 
presently make the following conclusions. 

1. We have spectroscopically confirmed, through the 
presence of 27 member galaxies, an over-density of 
galaxies at z = 1.7089. 

2. The 3.6/im photometric richness is measured to be 
A^gai=30±8. Applying standard, low-redshift, con¬ 
versions, the mass (within 500kpc) of the cluster 
is estimated to be ^ 3.8±1.2xlO^^M0. The ve¬ 
locity dispersion, calculated from all 27 members 
over within 1.8 Mpc of the BGG, is ddOj'j^Qg km/s, 
implying a lower mass of 8±3xl0^^ Mq. 

3. The Brightest Cluster Galaxy is coincident with an 
Ultra-luminous Infrared Galaxy with an estimated 
luminosity of Lir = 6.2±0.9xl0^^ Lq. The pres¬ 
ence of PAH lines in its infrared spectrum and its 
rest-frame NIR colours indicate this energy is dom¬ 
inated by star formation, with an AGN corrected 
rate of of ^ 860±130 Mq yr“^. This rate necessi¬ 
tates rapid gas deposition into the BCG system. 

4. Deep HST NIR imaging has revealed complex mor¬ 
phology within the core. A 66 kpc diffuse struc¬ 
ture, reminiscent of a tidal tail, appears to emanate 
from the BGG. This structure contains several 
(>10) clumps of objects along its length, aligned 
as “beads on a string”. The MIPS emission is co¬ 
incident with a set of red beads in the center of the 
tail. 

5. The overall arrangement of the system suggests 
that the star formation is being driven by gas-rich 
galaxy-galaxy interaction, possibly involving the 
BGG. If so, this would be a new example of BCG 
mass growth through wet-mergers at high redshift. 
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